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Executive Summary 

This report describes the analysis by MSE Technology Applications, Inc. (MSE) of the Electron 
Spiral Toroid (EST) concept being developed by Electron Power Systems Inc. (EPS). This 
analysis was conducted for the National Aeronautics and Space Administration s Langley 
Research Center (NASA-LaRC). The EST is described by EPS as a plasma shell of toroidal 
shape, in which electrons are in orbits around the toroidal axis. The current produced by the 
electron motion generates a toroidal magnetic field, which is confined within the torus shell. It is 
claimed that this plasma structure is very stable and can store vast amounts of energy, primarily 
as a magnetic field. The most detailed description is found in a theoretical study of the EST 
concept performed by the Massachusetts Institute of Technology (MIT) and published by EPS, 
herein listed as reference [I]. EPS also recently published the Phase I report for a NASA Institute 
of Advanced Concepts (NIAC) award, in which some of the specific claims about lifetime and 
energy density storage were made, denoted herein as reference [2]. The present analysis was 
originally based on these and another document describing the EST concept, which was also 
published by EPS [3]. However, dissemination of an early draft of this report led to a further 
exchange of information between MSE and EPS [4]; the release of two new documents [5,6]; and 
multiple revisions of the model by EPS. Reference [5] is actually a revised version of [3], and the 
second document [6] is an experimental report for the Defense Special Weapons Agency 
(DSWA). The analysis of this additional material is included in the final version of this report. 

Following the original description of the EST, the present analysis is mainly based on the cold- 
fluid equations for the plasma shell coupled to the self-generated magnetic field. The dynamical 
equations of the cold electron fluid, as well as the numerical solutions, are found to match those 
presented in [1], However, these equations only represent an incomplete subset of the plasma 
dynamics, since the ion motion was ignored. It is easily seen that the ion fluid is unstable, due to 
the repulsing electrostatic field required for electron confinement. A detailed analysis with the 
numerical procedure confirmed these results. The EST concept requires an excess of positive 
charges, and the potential well thus created is responsible for confining the electrons as well as 
the magnetic field stored within the torus. Although various scenarios for ion confinement were 
then suggested [4,8] by EPS and MIT, none of them is plausible. In fact, it is demonstrated in this 
report that this problem is fundamental and cannot be solved within the context of the proposed 
EST configuration. Several problems were also found with the results presented in [2] regarding 
the physical characteristics of a high-energy EST, which were found to be inconsistent with the 
scaling properties of the numerical solutions. The rates of energy losses through collisions and 
cyclotron emission were also examined and found to be in error by several orders ot magnitude. 
Furthermore, the analysis by EPS of the experimental data [6] does not provide proof of an EST, 
as defined in [1]. 

It is demonstrated in this report that the claims of absolute stability and large energy storage 
capacities of the EST concept have not been substantiated. However, there is undeniable 
experimental evidence of some type of plasma structures whose properties and characteristics still 
remain to be determined and which could potentially have applications tor space propulsion. 
However, more realistic theoretical models must first be developed to explain their existence and 
properties before applications of interest to NASA can be assessed and developed. 


I 


1. Introduction 


This report describes the analysis by MSE Technology Applications, Inc. (MSE) of the Electron 
Spiral Toroid (EST) concept being developed by Electron Power Systems Inc. (EPS). This 
analysis was conducted for the National Aeronautics and Space Administration’s Langley 
Research Center (NASA-LaRC). The company received a Phase I award from the NASA Institute 
for Advanced Concepts (NIAC) - grant number 07600-013. The information regarding the EST 
concept is somewhat sketchy, since there are no peer-reviewed publications. Our first analysis of 
the concept was completed in November 1999, and was based on the limited literature available 
at the time; primarily a theoretical description of the concept by C. Chen [1], the NIAC Phase I 
report [2], as well as an EPS document describing other aspects of the concept [3]. After the 
release of an earlier version of this report and a discussion with EPS investigators [4], further 
information was made available [5,6] and subsequently reviewed. Finally, a meeting took place 
with EPS and the Massachusetts Institute of Technology (MTI) [7] during which the latest, 
nonpublished information [8,9] was made available. Note also that references [3,5,8] are different 
versions ot a similar document; the revisions were made following the results of the present 
analysis. 

In the original ’ published” documents [1-8], it is claimed that a “stable” plasma, capable of 
storing vast amounts of energy, can be created without the need for external confinement 
(i.e., free ESTs). There is no exact definition of “stable”, but indications from the available 
literature [2,3] suggest lifetimes of the order of years. Since plasmas in laboratories are 
notoriously short-lived entities, this claim is quite spectacular. Other claims are similar in scope. 

For example, the EST specific energies could be of the order of 6xlO if> J/kg. By comparison, 
combustion ol I kg of hydrogen with oxygen would generate 1 .2x1 0 8 J, which is equivalent to the 
liberation of approximately 1.25 eV/atom. To reach the claimed level of specific energy, each 
atom would have to be able to store 650 MeV of energy (approximately 60% of the rest-mass of 
the proton itself) ! To date, only antimatter is known to have this level of energy density. The 
stored energy would primarily be in the form of a confined magnetic field, of the order of 
16,000 Tesla. This would lead to an internal pressure of I billion atmospheres, yet no external 
confinement method is mentioned. 

The plasma in the EST is basically contained within a toroidal shell, with the electrons at high 
velocity in the poloidal direction (i.e., in orbit around the centerline of the torus). The electron 
motion generates a magnetic field along the toroidal direction, which forms the greater circle of 
radius R r . A schematic ot the concept is shown in Figure I below (taken from [2]). 



Figure 1: Schematic of EST plasma. 
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The plasma is assumed to be strongly coupled. The analysis of plasma dynamics in [1] is based 
only on the cold-fluid model, which is certainly valid for strongly coupled plasmas. It is sufficient 
to say that thermal effects are assumed to be nonexistent or negligible, and therefore, they are not 
considered in [I] or in the present analysis. 


2. Fundamental Dynamics 

The toroid of inner radius r h and outer radius R r , shown in Figure 1, can be approximated by an 
infinitely long cylinder. The examination of the stability properties of the toroid can then be 
reduced to a two-dimensional problem. As described in [1], the plasma is considered to be within 
a toroidal shell, between the boundary radius r h and an inner boundary r, . The plasma is 
therefore contained in a cylindrical shell. This configuration is identical to the one described in 
[1]. Consider also the transformation from an orthogonal coordinate system ( x, y, z ) into a 
cylindrical coordinate system (r,<p, S), as shown in Figure 2. The associated transformation 
matrix R ’ describes a rotation by an angle (p , such that: 
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Figure 2: Schematic of coordinate transformation. The z-direction is 
out of plane. The toroid is approximated by an infinitely long cylinder. 


The equation of motion for a cold, incompressible plasma fluid is: 

^ + (u n • v)j„ = ^SL (E + U a x B) (4) 

at m a 


’ In the torus, the z direction becomes 6 , the toroidal direction. 
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where a = e,i is the index for the plasma component (electrons and ions). The dynamical 
Equation (4) can then be expressed in the cylindrical coordinate system, using the fact that the 
magnetic field has a component in the z direction only: 


B,u r +u r d r u r -— = —(£, + « B.) 
r m a 

d,u,„ +u r d r u m + - - - (Ey-u B 7 ) 


m 


(5-a) 


(5-b) 


a 


Since one is interested in steady equilibrium conditions, the time-derivatives and the radial 
velocity are assumed to be identically zero. The equilibrium configuration being symmetric, we 
also have E v = 0. The only velocity is in the q> direction, the only electric field is in the 


r direction only, and the magnetic field is in the z direction. Therefore, one can drop the spatial 
index for these quantities. Assuming a single species of ions of mass m l , we then have: 
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r 



r 


=—(E + u e B) 
m e 

— (E -i -UjB) 


( 6 ) 


The electric field is given by Poisson’s equation, i.e.: 

-d r (rE) = — (Z i n l (r)-n e (r)] (7) 

r £" 

At this stage, it is useful to introduce density scales , n ! and bring all the spatial dependence 
into normalized functions f e (r), /.(r) , such that: 


nA r )= n efe( r ) i Z i n i {r)=n l f i (r) and A/(r)=/,.(r) -f e (r) (8) 

A spatial scale is also introduced as the minor radius of the toroid, r h and a normalized radial 
variable defined as: 


p=— ; 0<p<l 


Introducing also the plasma frequencies co pe ,co {){ : 

i 

e~?i. 


<°~re = 


m e , 


c o 


, _ Z l e-n i 


"he,, 


Poisson’s equation becomes: 

eE 




(p)= 

m r 


jdp'p'Af (p') 


4>(p) 


The solutions of the equations of motion then become: 


“AP) = ~ 

«,(P) = 


ft) (r ± +4ft)^,ft>(p)/p 2 
± a /",v- 4 <^ > (p )/ p 2 


(9) 


(9) 


( 10 ) 


(H) 


where the definitions of the cyclotron frequencies have been used: 
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eB s ... o Zjtn, 

ft) = — ; ft),, = <5ft)„, with d = 

“ • m m, 


( 12 ) 


Finally, the magnetic field is given by Ampere’s law: 

d r B = -ejl n {n,f l u l -nJ'U e ) (13) 

Assuming that the plasma is nearly neutral, one can set ii, = n e , with any space-charge effects 
included in the distribution functions In that case, the evolution of the magnetic field 

becomes: 

c 0 : 


d,xo ce = (f i u i -/,«,) 


(14) 


Equations (I0),(1 1) and (14), combined with the definition (12), form a closed system and can 
thereby be solved. 


3. The EST Configuration 

According to ref. [I], the EST is characterized by an electron density function of the form: 


fAP) = 


\rj r if r, <r<r h 


(15) 


0 otherwise 

v 

The ions form a background density with the same spatial dependence as the electrons, albeit with 
a slight relative excess f , such that: 

In that case, 


and 


II 

'-K 

(16) 

0(p)= [ dp'p’Af(p') = — (/-l)(p-— ) 

r,lr h ^ ^ 

(17) 

•O St 

II 

1 

1 

^ 1-^ 

(18) 


Although there is no other explicit information regarding the ions in the EST model, it is assumed 
that they provide a fixed background since the ion equation of motion is not being considered in 
[1]. It should also be pointed out that there is a somewhat confusing statement in [2], where the 
parameters of a large-scale EST are listed (page 1 1 of ref. [2]). Although the total charge is 186 
coulombs, the system is called “charge-neutral”. Even though the relative deviations from 
neutrality would appear small , the densities considered in these intense EST configurations lead 
to a large amount of electrostatic energy. In fact, it is precisely this electrostatic charge that is 
responsible for containment of the plasma and storage of the intense (16,000 Tesla!) magnetic 
field. At sufficiently large distances from the EST, the electric field given by the EST can be 
approximated by the zero th -order term in the multipolar expansion: 

Q 


E «- 


( 19 ) 


AtieR' 


For example, the last theoretical case considered in [ 1 ] (page 16 ) uses a deviation /-I - 5 x 10 


-7 


5 


which (for the 186 Coulombs quoted in [2]) gives a field of £ = 1.7xl() 12 Vim at 1 m distance. It 
is not known how the developers of the EST intended to prevent electrical discharges from 
spontaneously occurring between the EST and its surroundings, given this level of field intensity. 


Nevertheless, the theoretical analysis of the EST concept, as described in [1], can be continued. It 
is useful to introduce further normalization parameters, such that: 


and 


v e=- 


li e 

°>pe r \ 


^ pe^\ 


(20-a) 


Cl=^ ; Cl =^~ 

n‘ ^ ’ Cl 


CO n 


(Q„ 


(20-b) 


Instead of normalizing the distance to the boundary of the toroid (Equation (9)), one can scale the 
distance with respect to the inner radius of the electron shell, as in [1]: p = r/r t . In that case, the 
electron velocity is given by: 


£2 + 

ce — 




Ll;. e+ 4(f-\)j(\-±) 


( 21 ) 


Introducing the coupling parameter: 


a = 




and since the ions are fixed, the equation for the magnetic field can simply be written as: 


dCl 


dp 


— Gr- 


and therefore: 


dCl Ct , _ a 

dp 


Q„±/q,;+4(/-1)-E(1-±) 


(22) 


(23) 


(24) 


Equations (21) and (24) describe the evolution of the system and are completely identical to 
Equations (23) and (24) or ref. f I ]. 


Therefore, the dynamical system of the EST (as described in [1]) has been recovered. This is a 
nonlinear system, which requires a numerical integration procedure. Therefore, a simple 
numerical procedure tor the solution of the system (21, 24) was written, and numerical solutions 
were obtained for the same parameters as those given in [1]. For example, the normalized 

solution for a = 1.6, / = 1.02, and r h !r x =1.01 is shown in Figure 3. Note that the boundary 
radius r h used here is equivalent to the r h2 radius used in [1]. The solution presented in Figure 3 
is completely identical to the corresponding Figure 2 of ref. [1]. All the cases computed in [1] are 
listed in Table 1 and have been recomputed with our numerical integration code. These solutions 
are shown in Figures 3 through 8 and are identical with the results presented in [1]. However, 
there are no figures corresponding to the last case (#4) in [1]. Special attention has been given to 
the solution accuracy, and the current numerical scheme uses a stretched grid near the edges of 
the plasma shell, such that a smooth solution is guaranteed in that region. 

Using the parameters listed in Table 1, the physical solutions can be obtained from the 
normalized numerical solutions. Ot special interest are the peak electron velocity (at the outer 
edge of the toroidal shell), the peak magnetic field (at the inner edge of the shell), and the total 
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number of electrons (and ions) in the shell. The latter can be approximately obtained from the 
electron density and the volume of the toroidal shell, neglecting the variation of the density 
(n e r [ /r) within the thin shell (r M -r,«r, ). Since r, and a are given, the plasma frequency can 
be obtained from Equation (22) and the electron density from the definition of the plasma 
frequency (9). The volume of a thin torus ( R T »r h ) is approximately given by: 


Volj ~ 2n 2 R T r,, 2 

and since the plasma is contained wtthm the shell delimited by r t ,r h2 , its volume is: 

Vol p =2 Jt 2 R T (r; 2 -r{) 

And therefore, we approximately have: 


N , = Vol, 


( 9 V 

ac~ m ' 


n 


V 


(25) 


(26) 


(27) 


r{ e~ 

A more accurate value of the total number of electrons A/^can be obtained by integrating the 
density profile during the numerical integration. This is accomplished in MSE’s code. 


The results for all cases considered by Chen in [I] are shown in Table 2. Note that the velocities 
and magnetic fields can be either positive or negative, depending on the sign of the solution taken 
in Equations (21) and (24). Only the absolute magnitudes are listed in Table 2. With a few 
exceptions, the current results compare very well with Chen’s. The discrepancy in N t for the first 
two cases is obviously an accounting error in [1]. Our results are consistent with the 
approximation (27). Another discrepancy is in the field value for Case #4. Since the normalized 
solution was not shown in [1], it is difficult to assess whether the error lies in Chen s numerical 
solution or in the rescaling. Since the electron velocities agree very well for this case and the 
numerical accuracy is better for the field than for the velocity profile (the velocity gradient 
becomes very large near the boundary at large a, thus necessitating a highly stretched grid), one 
can reasonably assume that Chen’s numerical solution is also accurate. Therefore* it is assumed 
that there is a simple scaling error in [1], As will be discussed below, our values of the magnetic 
field are seen to obey strict scaling laws, thus adding validity to the current results. 


Case # 

a 

/-I 

rjJr, 

r, (m) 

Rj (m) 

i 

1.6 

2x1 O' 2 

1.01 

4.95x1 O' 4 

2.5X10' 3 

2 

10' 

1 x 1 O' 4 

1.01 

1.00x1 0' 2 

0.10 

3 

2xl0 4 

4x10 s 

1.05 

1.50x10 2 

0.15 

4 

4x10 s 

5x1 O' 7 

1.05 

1.50x10 2 

0.15 


Table 1: List of cases studied by Chen in [1]. 


Case 

V' (r 2 ) 

(m/s) 

Be(f, 

(T) 


r. 



1^# Smmt $ j| | 




current 

i 

5.37x10 s 

5.39x10 s 

6.60x1 O' 4 

6.64x1 0" 4 

2.24x10'° 

4.50x1 0 10 ’ 

2 

9.49x10 s 

9.48x10 s 

1 .49x1 0" 2 

1.48x1 O' 2 

5.6lxl0 14 

1.1 1 x 1 0 1 5 7 

3 

6.00x1 0 7 

5.99x1 0 7 

0.69 

0.70 

1.68xl0 17 

1.67X10 17 

4 

9.5 IxlO 7 

9.48x1 0 7 

21.2 

15.58 * 

3.36XI0 19 

3.34xl0 19 


Table 2: Comparison of our numerical results with those of Chen in [ 1] . 


: it would appear that the total number of particles in 1 1 1 have been erroneously divided by 2 in these two cases. The 
correct number can easily be obtained from the torus volume and the electron plasma frequency, using Equation (22). 
? Notice the discrepancy in magnetic field in this case, while the velocity is in good agreement. 
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1.04 
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0.00 L - 
1.00 


r/r i 


Figure 6: Normalized cyclotron frequency for Case #3 (compare with Figure 5 of [ I ]). 
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IL/cu, 



Figure 7: Normalized velocity for Case #4 of Table 1 (symbols are 
plotted every 5 grid points). 



Figure 8: Normalized cyclotron frequency for Case #4 of Table I (symbols 
plotted for every 5 grid points). 
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4. Stability Problems 


Both the dynamical equations and the numerical solutions of ref. [I] have been recovered, which 
would seem to bring some validity to the EST concept. However, it is crucial to examine the 
dynamical behavior of the complete system of equations including the ions'. A key assumption 
barely mentioned in [1] is that the ions form a fixed background. However, since the magnetic 
field vanishes at the boundary of the toroid, there is no confining force being applied to the ion 
fluid near the surface. Since the toroid is also positively charged (the electrostatic force provides 
the confining force for the electrons), there is a strong repulsion pushing the ions outward. 
Therefore, there is no mechanism for keeping the ions inside the toroid. They would very rapidly 
“evaporate”, thereby leaving the toroid absolutely charge-neutral. This process would also 
eliminate the confining force for the electrons, and the whole toroid would then rapidly expand, 
unable to contain the magnetic pressure. To look at the ion dynamics, consider Equations (5-a,b). 

For no initial azimuthal (i.e., poloidal) velocity (u^sO), and since E r >0 (/> 1 for electron 
confinement). Equation (5-a) clearly yields d,u r > 0, and the ion gas is rapidly expanding. The 
fixed ion background assumption is therefore completely invalid for ion-motion time scales. 


One could then argue that a stable solution could still be achieved, but with the ions in motion. It 
is clear that as ions gain a radial velocity, the Lorentz force will also act to provide them with an 
azimuthal velocity (see Equation (5-b». In that case, one must examine the stable solution of the 
ion fluid with an azimuthal velocity. The ion solution has already been described in Equation 
(II). Using the same normalization (20-a,b), one obtains (for the same density distribution): 


Vi 2L 


p [- 4 , ± 

Using ( 1 2), this can also be written as: 

P 

2 


v, ± 1 /ni-*</-i)£0-£> 


(28) 


(29) 


A stable solution is possible if the discriminant remains positive, i.e., if: 

A i = Cl;,-'¥ 2 >0 (30) 

with 



Obviously, since the field vanishes at r = r h2 , Q ( , =0 and the condition cannot be satisfied. In 
fact, there will be a region of finite extent near the outer boundary, where this condition cannot be 
satisfied. The size of this region depends on the parameters a and /-I . Figures 9 and 10 show 
plots of both Cl;, and X P 2 for Cases #3 and #4 of Table 1 . The allowed region where a physical 
solution can be found for the ions is determined by the distance between the origin ( r = r, ) and 
the point at which the two curves cross each other. For Case #3 (Figure 9), it appears that most of 
the plasma in the shell would be unstable to ion “evaporation”. For Case #4, this region is much 
smaller. One could argue that for the right choice of parameters, the unstable region would be 
reduced to a “skin”. However, this does not solve the problem since the ions would always be 

emitted from the surface of the EST. As the relative charge ratio f drops, the electrostatic 
confining force would also drop, and the torus would rapidly expand under the pressure ot its 
own magnetic field. As the EST surface increases during expansion, so would the rate of ion loss. 


1 1 


0.016 



Figure 9: Comparison between the two terms of Equation (30) determining 
extent of physical solution for the ion velocity. Results for Case #3 of Table I . 



Figure 10: Comparison between the two terms of Equation (30) determining 
extent of physical solution for the ion velocity. Results for Case #4 of Table 1 . 
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The process is self-accelerating and would eventually lead to the explosive liberation of the stored 
energy. However, there are still further problems that prevent a stable solution. As the A, term 
becomes positive and a physical solution for the ions can be found, the dynamics are no longer 
governed by Equations (21) and (24). In fact, the full expression for the current must be used 
(including the ion velocity) and the evolution of the magnetic field is now given by Equation 
(14), or its normalized equivalent. Solving for the complete system, it is found that no physical 
solution can be obtained. Detailed examination of the results reveals that the problem has nothing 
to do with numerical errors. It is first necessary to point out an important aspect of the dynamics 
of the problem. There are two possible solutions for the electron velocity, depending on the sign 

chosen in Equation (21). If the “+” sign is chosen (and t>+>0), the field takes on negative values 

(t2*,<0). If the sign is chosen (and i+<0), the field is positive (Q^>0). This can easily be 
seen on Figures 3 and 4. In either case, the dynamics are such that the absolute value of the 
electron velocity is given by a difference between two terms: Cl ie and the square root term of 
Equation (21). Therefore, using the definition (31), one can write: 




I-.1 + 


<5 4 *- 

al 


(32) 


However, at the crossover point where the ion fluid can have a physical solution, Q (( , ~ y , and 


therefore: 


(33) 


in that region. 


On the other hand, the ion velocity is approximately: 

v,—s£a„ « 4 > 

since the square-root term in Equation (29) vanishes at the crossover point. Therefore, at the 
crossover, the ion velocity is approximately twice the magnitude of the electron velocity and with 
the opposite sign. The rate of evolution of the magnetic field immediately changes sign, and the 
magnitude of the field starts decreasing as one moves towards the inner boundary. This also 
implies that the discriminant A, given in Equation (30) immediately ceases to be positive, and the 
ion fluid has no physical solution. This behavior is exactly confirmed by a careful examination of 
the numerical integration procedure as one attempts to find a solution that includes the ion fluid. 
Therefore, even if the parameters are chosen to concentrate the dynamics near the outer edge of 
the shell, the complete dynamics reveal that no stable solution can be found. 

It is also useful to estimate the field required to enforce a stable solution for the ion fluid. This 
could be estimated by requiring that A, = Oat the boundary r = r h2 . For example, taking Case #3 
of Table 1, one finds that there must be a magnetic field at the outer boundary r b2 ~ 1.01 r, on the 
order of: B ex , =0.9 T . This is comparable to the field stored inside the EST for this case. 
Therefore, even if one were to provide an external confinement field, its value would need to be 
of the same order as the one stored inside the toroid. 


As pointed out by Chen [4], the discussion of ion stability can be avoided altogether by invoking 
the Virial theorem ([10], pp. 72-74), which states that a plasma cannot be self-confined. However, 
we believe that the examination of the dynamical equations elucidates the origins of the self- 
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confinement problem better than maybe simply quoting a theorem. At the same time, our analysis 
indicates that an externally confining field would have to be of the same order as the stored field. 
Although this may appear obvious from energetic considerations, the result is important because 
it can be generalized to any method of confinement including external mechanical pressure. The 
magnetic energy stored inside the torus is in direct relation with the electrostatic potential energy 
generated by the excess charge ot the ions used to confine the electrons. In turn, to confine the 
ions would require an energy density (i.e., a pressure) of at least the same order of magnitude. 

In the introduction of ref. [1], it is mentioned that the model is valid only for time scales short in 
comparison with the ion motion time scale (e.g. inverse ion plasma frequency); the problem of 
ion confinement is therefore not addressed. This is the only mention of this important time-scale 
limitation in all the provided documentation. On the other hand, several EPS documents [2, 3,5,8] 
state that the physics of the EST were derived and verified, the EST was stable with no external 
field required for confinement, and there was “no known or obvious normal occurrence that will 
lead to instability”. There can be no ambiguity in the definition of “stable plasma” since the 
described applications require lifetimes of many hours, certainly much greater than the ion 
motion time scale. For example, references [3,5] explicitly mention an EST lifetime of I0 6 hours 
[3,5] in a vacuum environment. 

It is then clear that the free EST, as described in [1], cannot achieve long-term stability. Further 
discussions with Chen [4] confirmed that mechanical pressure of some sort is required for ion 
confinement. The realization of a high-energy EST in a vacuum environment is therefore 
completely unfeasible. In a follow-on discussion [4], EPS officials implied that the references 
were actually for a “partial vacuum” only. This is not a valid argument, since the term “vacuum” 
can be used to describe any environment where the pressure is much lower than atmospheric and 
therefore, negligible compared to the internal pressure (gas-dynamic or electromagnetic) of the 
high-energy EST. At any rate, it is extremely difficult to imagine how high-energy ESTs can be 
confined by a mechanical pressure many orders of magnitude smaller than the stored magnetic 
energy. The example described in [3,5] mentions an EST with 100 MJ of stored energy, 
corresponding to an energy density of approximately I0 M J/m 3 (i.e., a magnetic pressure of 
2 million atmospheres). The example in [2] corresponds to internal magnetic pressures of 
I billion atmospheres, yet one should expect these high-energy EST configurations could be 
absolutely stable in near- vacuum conditions! To find a mechanism that can achieve this 
extraordinary feat would require dramatic changes in the fundamental laws of physics. The 
problem of ion confinement will be further discussed in Section 7. 

5. The EPS Theoretical Model 

Before discussing further the implications of the cold-fluid model, it is worth investigating 
another theoretical model of the EST proposed by EPS, as described in [5,8], This is an attempt a* 
constructing a "discrete” model using a summation over a finite number of particles. The ions and 
electrons reside in separate shells (the ions form the inner shell), composed of a large number of 
loops of radius r 0 and a large number of particles along each loop. The spacing between 
electrons on a loop is d e , and therefore, the number of electrons in a loop is: 

M ~ nr o 

^ {e\loop) ~ j (35) 

The number of loops in the toroid is approximately ( R T » r () ): 
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(36) 


N 


(loop) 


2nR T 

(k 0 cl e ) 


where k 0 d e is the spacing between loops, with k„ ~ 1 . The total number of electrons is then: 

4n 2 f 0 R T 


N e 81 (elloop)" (loop) 


(37) 


k ( ,d; 


Note that the EST surface area is also given by: 

Area = 4n 2 r„R T . ^8) 

The EPS investigators first compute the force acting on an electron in a single loop. Since the 
distance between electrons is given by the chord d e = 2r n sin(50 12) ~ r ft 86 , where the angular 
distance is of course 86 = 2n / N (eVoop) , the radial compon ent of the electrostatic repulsion force 

from all electrons in a loop is then: 


F. = 


\KE, 


f l/wyi ) 

■2 t 
// = 1 


/ 2 


sin(/i<50 / 2) 


dl 


(39) 


where d n = 2 r„ sin (n86 / 2) . The EPS investigators then use the approximation of small angles to 
write the final answer as: 


Fe 


4ne 0 r 0 d e 


,1 1 1 . 

1 + — + - + — + 
2 3 4 


(40) 


Strictly speaking, this is a diverging series and a potential problem for large i oof)) . EPS 

estimates the summation to be «5, a rather arbitrary choice. Using this small finite value for this 
summation, EPS then expresses the force of repulsion due to all other electrons on the loop as: 


F. =C' 


Any r d 


(41) 


where C is a constant. Note that this is the force acting on a single electron, and that this result is 
obtained in the limit of large N {eVoop) . EPS then uses this expression in the balance of forces. 

However, this is a serious error for several reasons. 

First, the expression for the force itself (41 ) is wrong, because the truncation of the summation to 
small numbers is invalid. The summation in (40) is logarithmically divergent, and the expression 
for the force should then include a term log(r„ / d e ) , which is not negligible. 


Second, and even more importantly, this expression does not take into account the contribution 
from neighboring loops. Since the distance between loops is assumed by EPS to be of the same 
order as the distance between electrons in a loop ( k„d e , with k a ~ 1 ), this contribution is far 
from negligible and is, in fact, of the same order. The correct evaluation of the force is shown 
later in this section, but an estimation of this contribution can be given below. 

To include the electrons in neighboring shells, one can replace each electron on the loop by a tube 
extending in the z direction (the major radius of the toroid). This tube has a diameter d e , and the 

effective density of electrons in this tube is (n 'd;k 0 d e / 4) _l . The solution of Poisson s equation 
for this charge configuration yields the following field strength at a distance d from the tube: 

E = !__ (42) 

ne 0 k„d e d' 
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Using this value for the radial force on an electron and summing over all tubes along the loop: 


F. = 


ne (1 k 0 d e 


8 ^g^ sin (nS0/2) 


(43) 


n=\ 


This time, the sum can be evaluated since the sin(/i<50 / 2) cancel out. The summation becomes: 

1 N, 




n=l 


2 r. 2r„ 


“ (elltiop) _ K 

_ 4Z 


(44) 


The force on the electron then behaves as d e 2 , not as ( r ( ,d e ) 1 and is therefore much larger. 

The errors made by EPS are then further compounded when computing the force due to the 
magnetic field. This is expressed in [5,8] as: 

evpt'JN 


F„, = ev x B = - 

2nR, 


(45) 


where i = ev c / d e is the current in a single loop, and R T is the radius of the toroid. There is some 
confusion as to the nature of this field. In [5,8], the magnetic field is taken as an initiating field 
caused by the arc current used to create the EST. In that case, it is not clear why the arc current 
should also be the same as the loop current, unless all the arc current goes into the toroid. In that 
case, the magnetic field created by the current loop should be: 


B = 


2/rr 


(46) 


and the toroid radius Rj should not appear in these expressions. There is also confusion about 
the meaning of the number N appearing in (45). This should actually be the number of loops: 

2nR r 


N 


loops 


k„d. 


(47) 


but in fact, EPS uses the product N = N loop N (eUoop) in Equation (45). The resulting expression in 

[5,8] yields the force acting on all electrons in a loop instead of on each electron. The 
comparison with the electrostatic forces (Equations (41) or (43)) is therefore completely invalid. 


The proper way to evaluate these forces can now be described in the same limit of large N 
Assume that the electrons form a cylindrical shell of thicknessdr . Integrating Poisson’s equation 
(7) for the electron component only, the electric field acting on each electron becomes: 

£7 en e , 

E = dr (48) 

£<> 

where n e is of course the electron density in the shell. Using the total number of electrons (37) 
and the shell volume (27rr„ )dr(2nRj ) , one can easily arrive at the final expression: 


£ = € — and F e =-^ 

£ 0 k ( ,d e 


(49) 


Notice that this result is independent of the assumed thickness. The field due to the ion shell can 
be computed in a similar fashion and added (linear superposition) to yield the total electrostatic 
force on an electron. The magnetic field inside the toroid is approximated by a solenoidal field: 


... 


B = U„i 


dx 


■ AV 


k () d e 




16 



or: 


B = — ~r=v e 


k 0 d; 


e„k„d: 


Note that B, =-B . The magnetic force on each electron is then: 


- 


k ( ,d] 


i j 2 2 


(51) 


(52) 


Note that the ratio of magnetic to electrostatic forces is v~ / c ~ , as expected. These expressions 
can now be used to evaluate the balance of forces, and the contribution of the ions can be seen to 

be / times the force due to the electrons (49). The balance of forces for the electrons is: 

(f-\)-^ = ev e B + m e ^ (53) 

r o 


Another application consists of examining the stability of the ion shell inside the electron toroid. 
By symmetry, there is no electronic contribution to the electrostatic field, and the balance of 
forces (F lon +F rot =F mag ) on each ion (if possible) reduces to the following equality: 


•» 2 
- v 

f + rn, — = eVj B 


where d , ~ d e is the ion spacing. Let us now introduce the ion Larmor radius: 


(54) 


r L{i) =- 


m , v , 

eB 


‘ £ d 2 ^ 

e 


rw,c\) 


( V ^ 


(55) 


where we have again used Equation (51) for the magnetic field. The centrifugal pseudo-force in 
Equation (54) can now be brought into the form: 


m, 



S 

Cl 

( r ) 

r LU) 

v,v e 

r 

y£ () k () ct^ J 

r 

\ J 

c~ 


The balance of forces for the ions (54) then implies: 


(56) 


/ 


£ 0 k 0 d~ 


V 


^ £ ( tk 0 d e ^ 


’Hi) 


\V;V e 


(57) 


Both the centrifugal and magnetic forces are of the order of VjV e /c and can be neglected 
compared to the electrostatic force. Therefore, there is no stable physical solution for the ion 
shell. Of course, this was also demonstrated in Section 4. 


It is also interesting to examine the consequences on a mechanical confinement scheme. Suppose 
that the ion stability is provided by a mechanical (gas or containment vessel) pressure. In that 
case, the confining force applied to each ion must balance the electrostatic repulsion. The balance 
of forces then becomes: 

2 v 2 

f “ + Wlj — = eVjB + F(conf\i<m) (58) 

£ <>k()d j i 


Neglecting the terms of order v,v e /c 2 , one can then express the required confinement force for 
the entire EST: 
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(59) 


P ( conf I EST) p ( cirnf lion) ^ {ion ) 


Y 


y £ (>k,,d j 


f 


4 n\Rr_^ 

d] 


where the total number of ions is / times the total number of electrons, given by Equation (37). 
The total confinement force can be divided by the area (38) to give the confinement pressure: 


P amf = f 


£ 0 kldjdf 


(60) 


This expression can be manipulated further. Using (51), the electron spacing can be expressed as 
a function of the stored magnetic field: 


; ^2 eV c 

M, =/*„- r- 

D 


In the approximation of nearly identical spacing d; = d e , we therefore have: 




£„k~d-d; ~{e o n o )e 2 -t- 


t B 2 ' 


-i 


2/U, 


The confinement pressure becomes: 


P con,=f 


2 £- 


U B 2 \ 




We can also express the ratio of speeds as a ratio of kinetic to rest mass of the electron: 


2c~ m e c~ 

and the final expression tor the confining pressure becomes: 

511 


P amf “ f P ma K 


E e [keV ] 


(61) 

(62) 

(63) 

(64) 

(65) 


Note that f — 1 and unless the electrons are ultrarelativistic (in which case this analysis must be 
modified), the required confinement pressure needs to be several times larger than the magnetic 
field pressure inside the torus. Therefore, even if one could envision applying a confinement force 
on the ions , the required magnitude of this force is such that the concept of energy storage (i.e., a 
large internal magnetic field) is completely impractical. The relationship (65) is very useful 
because of its generality and independence with respect to the EST dimensions. 


We conclude that the theoretical model proposed by EPS is severely flawed in many aspects and 
does not provide a stable solution. The only correct theoretical model at this stage is the cold- 
fluid model, which is described in [I] and also in the present analysis. However, Chen’s model is 
incomplete, since the ion motion has not been considered. More complex descriptions of the 
plasma suggested by Chen (e.g. phase-space approach) would not change the fundamental 
properties of the model, and would therefore not solve the problem of ion confinement. For that, a 
completely different physical model may be necessary. 


6. Scaling Properties 

Although it has been determined that the EST (as described anywhere in the literature [1-8]) is 
not a stable structure, there are other claims to be evaluated. It is first necessary to examine the 


It is not clear how BPS proposes to apply a mechanical (e g. containment vessel or gas pressure) force to the inner shell of the ions 
without a fleeting the outer shell ol the electrons. As is shown in Section 7, it is not possible to prevent collisions between electrons 
and other particles if put into contact. 
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scaling properties of the structure, which will help in the examination of the claims of [2]. The 
incomplete solutions, found with only the electron dynamics, are analyzed as function of the 

“structural” parameters oc and f—\ . The dimensions of a large-scale EST are chosen as follows: 

= 0.15 m, r h2 /r y = 1.05, R T = 1.5 m 

It is therefore a scaled-up version of Cases #3 and #4 of Table 1. A series of numerical solutions 
can then be obtained for various values of, Of and /—I and characterized as a function of the 
maximum magnetic field 5 niax (at the inner boundary r, ) and maximum electron kinetic energy 
E e max (at the outer boundary r /j2 ). The results are shown in Figure 1 1 for a constant maximum 

field, chosen to be approximately 16,000 Tesla. The results for a constant maximum electron 
energy (chosen to be 12.5 MeV (sic)) are shown in Figure 12. 

f-1 



Figure 11: Scaling relationship between coupling parameter a and 
charge parameter /-I for constant magnetic field. 



Figure 12: Scaling relationship for constant electron energy. 
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One should point out that this choice of electron energy is completely arbitrary and is only 
intended to demonstrate the obvious scaling relationship shown in Figure 12. A large electron 

energy allows one to cover a wide range of values of the a and / parameters. Another choice, 
such as the 10 keV used below, would give a line parallel to the one shown in Figure 12. The 
physical answers for the magnetic Field and electron energy would also depend on another 
“structural” parameter, the relative shell thickness (i.e., the ratio r h2 l r x ). For fixed values of the 

coupling and charge parameters (a = 10", /- 1 = 10 8 ), the maximum field and energy are 
computed for several values of the ratio r h2 /r t , and the results are plotted in Figures 13 and 14. 


a=10 11 , f -1=1 O' 8 



Figure 13: Scaling of electron energy versus shell thickness. 


cxrrIO 11 , M=10' 8 



Figure 14: Scaling of magnetic field with size ratio. Both data and 
approximate fit shown. 
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It appears then that the electron energy scales exactly with the normalized shell thickness 
r h2 lr -\ . The scaling behavior for the magnetic field is more complex, and the formula shown in 

Figure 14 is an approximation. Since the peak electron energy is obtained at the outer edge 
(where the field vanishes), a simple scaling relationship can be obtained for the velocity. From 
(22), we have at the boundary: 


V „ <*= 


w 


( r, 


0/2, 


bl 

r \ 


1--^- 

n>2 


converting to physical variables and expressing the energy ( E 


0/2 


_ i 


( 66 ) 


, e -Em e v ~ ), we finally obtain: 


E 


t'.max 



(67) 


which is confirmed by both Figures 12 and 13. Combining all the data, the following approximate 
scaling relationships can be obtained: 


^v.nux ~ 2.56x10 


3/ 


•l> 2 

V 6 


a 
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\l/2 


keV 


Tesla 


( 68 ) 

(69) 


Although the electron energies in the EST can be quite high, it is assumed that the temperature is 
very low (i.e., the electrons are monoenergetic). In that case, the coupling parameter can be large: 


r = — >170 , with R 1 (70) 

kT Ane„R 

This implies that T ~ 2“ K and presents another engineering challenge. For example, refs. [2,3] 
describe EST formation with an “injected” electron beam. It is not clear how EPS proposes to 
generate an electron beam and capture it in a shell with a relative energy spread of less than 10 . 
There is also some confusion regarding the EST production method; in [5] an electric arc is 
mentioned instead. The arc should have an even higher energy spread. 


7. Energy Storage and Lifetime 

The EST parameters listed in [2] are different from the cases studied above and in [I], and since 
there is some missing information, one cannot extract the exact relevant quantities. Nevertheless, 
the right orders of magnitude can be obtained. The relevant parameters are listed in Table 3. 


Diameter, major, r T : 

0.63 m 

Diameter, minor, r„: 

0.21 m 

Electron energy: 

lOkeV 

Total Mass: 

1.86x10 6 kg 

Magnetic Field: 

16,000 Tesla 

Energy Loss, Radiation: 

1.16 J/sec 

Energy Loss, Collisions: 

8.9x1 0 ,0 J/sec 

Total Energy: 

1.9xl0 12 J 


Table 3: EST parameters listed in [2], 
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It is not clear how these values have been obtained, and there are some inconsistencies. First, the 
list mentions diameters, while symbols for radius are used. Assuming a hydrogen plasma, the 

total mass gives a total number of ions (and electrons): /V, = N t , ~ 1 . 1 2 x 1 0 21 . Since there is no 
information concerning the shell thickness, it is necessary to assume a value. Choosing 
r h2 = 1.05r,, the plasma volume is of the order of 6.7xl0' ? m, and the electron density is 

therefore: « t . ~ 1 .7 x 10‘ 3 m \ The plasma frequency is then (o pe = 2.3x1 0'\y~‘ , and the 


corresponding coupling constant (assuming the inner radius is r,=0.2l m) must be a ~ 2.6x1 0 8 . 
The plasma frequency also corresponds to an equivalent field of m e (Q pi ,le~ 132 T. This is in 
marked contrast with the claimed value ot 16,000 T given in [2]. Although some numerical 
solutions (not shown) have a peak value of the normalized cyclotron frequency £2 >1 , they have 

all been found to be of the same order. This is where the scaling relations can be useful. Inverting 
(68,69), one finds: 


[j±\ 

I O' 8 


= 0.426 


gnu g 

AT) 


In (r 2 /q ) 


r 2 /r { - 1 


( a 8 

, 10 " , 


= 2.155xlO -5 




[ln(r 2 //| )] 2 


(71) 

(72) 


Considering two typical values for the relative shell thickness, 1.01 and 1.05, and using the 
claimed values of E m% and £ niax in Table 3, one finds: 


rjr, 

a 

/■i 

1.01 

5.57x1 0 15 

7.0xl0 16 

1.05 

1 . 16xl0 15 

6.7x10 16 

Table 4: 1 

typical structural parameters 


satisfying energetic values in Table 3. 


A search in parameter space in the cas er 2 /r, =1.05 yielded the “best” numerical solution (i.e., 

one with £, . niax =I0keV and fi nHX = 1 6,000 T ) for a -1 .9x1 0 15 and /-l=5xl0" 16 , very close to 

the values predicted by the scaling laws. These values for the coupling constant are 8 orders of 
magnitude greater than the one derived from the mass of the system! Furthermore, these values 

correspond to an electron density n v >6xl0 21 ' m~' . This is higher than solid-state densities (the 

number density for solid copper is 8x1 0' 8 , while it is 3x1 0' 8 for water)! Any error in the scaling 
relations (62,63) would reside in the scaling of the magnetic field with respect to the shell size. 
As seen in Table 4, the effect of size is limited. One must then conclude that the values of 
magnetic field, energy, and mass given in [2] are incompatible with each other, as well as with 
the numerical solutions . 

The next step in the analysis is to examine the rate of collisions in the system. Since the electrons 
and ions are not separated in the model described in [I], the electrons will loose energy through 
Coulomb collisions with the ions. The only way these collisions could be prevented is if the ions 
were moving at the same velocity, in which case there would be no current generated, as well as 
no magnetic field. An expression for the rate of energy loss is found in [II], p.4 1 : 


It is emphasized here that the problem is not so much with the absolute scale of the density, although astonishing in its 
own right, but with the severe inconsistency of the parameters. 
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(73) 


dE„ 


Z}e A b, In A 


= —n, 


H(x,m i lm l ,) 


dt ' 4 

where E e is the energy of the electron, b i = (m, / 2kT t ]) u ~ , and In A is the Coulomb logarithm. 
For electron-ion collisions: 

In A = 31.3--jln(n 1 ) + ln(7]) (74) 

The function H is defined by: 

H(x,p) = erf{x)lx-2n~ ul (\ + /t)exp(-x 2 ) (75) 


with 


X = b,u t , - 


f n h E^ 

mj, 


1/2 


(76) 


with E e and T, expressed in the same units. For the parameters listed in Table 3 and assuming a 
hydrogen plasma at T, ~ 2" K , we have x * 3.4x 1 0 s . For large values of x, H = 1 / x . Using the 
low value of density derived from the mass (n, ~ 1.7xl0 23 ) and In A = 7.35, the final result is: 

^ * 10 12 J/s ~ 7000 keV/s (77) 

dt 

This is the rate of energy loss for each electron. To find out the rate of energy loss for the entire 
torus, it is necessary to multiply by the total number of electrons, which is approximately 

N e = 1 . 1 x 10 21 . Thus, the total rate of energy loss is approximately: 


> 10 9 J/s for the entire EST 


(78) 


Expression (77) implies that it takes approximately 1.4 msec for each electron to loose its energy. 
This corresponds to 300 orbits around the minor circle. Although the mean-free path is large 
compared to the dimensions of the system and the plasma can be considered “collisionless”, the 
neglect of collisions is valid on the time-scales characteristic of the plasma dynamics. However, if 
one looks at long-term stability, the absolute rate of collisions must be considered. Another way 
to evaluate the rate of collisions is to look at the Rutherford scattering cross-section, i.e., the 
single binary collision for a deflection angle of 90° or more. For example, using expression 
(6.1.15) of ref. [12] or (9-10) of ref. [10], the cross-section for such an event is approximately 

a,*, = 10 26 m 2 . This leads to a frequency of events 10 5 s' 1 Since only a fraction (=28) of the 
electron energy is transferred during each such collision, the time required for complete 
relaxation of the electron energy becomes 10 2 sec. However, at these plasma conditions, the 
impact parameter for 90° deflections p ^ is very small. The ratio of large to small angle collisions 

is then [10]: 


P 90 


V 


10“ 7 «1 


(79) 


and the expression (73) for multiple collisions is more appropriate. Therefore, one expects the 
rate of energy loss to be given by Equation (78), but it could be much larger if the density were to 
be consistent with the scaling relationships obtained from the numerical solutions. Note that EPS 
claims a rate of I O' 9 J/s for the entire EST, a difference of 18 orders of magnitude'. 


To emphasize the consequences of this, consider the following. Assuming the rate of energy loss 
is indeed given by the value in Table 1, the rate of energy loss per electron becomes: 

^4- = 8xl0“ 31 J/s (80) 

dt 

for the entire system. 


23 


This rate would be obtained for a number density of the order of 10 23 m 3 . Since the rate scales 
linearly with the number density of ions and if the number density was of the same order as the 
number density of copper , the rate of energy loss per electron would be: 

dE 

— ^ = 3.8xl(T 25 J/s = 2.3 x I O' 9 keV/s (81) 

dt ' ' 

It would therefore take 4xl0 4 sec for the 10 keV electron to loose its energy, which implies that 

such an electron beam would be able to go through 4xl0 13 km of solid copper without stopping! 

One could argue that the loss-rate by collisions (given by Equation (73)) is a local quantity and 
should be evaluated through volume integral, rather than using global and maximum values for 
the density and electron energy. Although the density varies as r,/r, the range of variation is 
small (since r^/r 1.05 ), thereby providing at most a 10% change. However, the variation of the 

electron energy is much more rapid and could make a difference. Therefore, the numerical 
integration procedure was modified to include a volume integration of the collisional losses. One 
should point out first that for 2000 , the function H( x) has a maximum near x~3 and can 

be approximated as Vx beyond that. Using local values (jc> 3) only in the integration, one finds 
that the volume integration gives a result within 25% of the result (78) obtained from global and 
peak values. Using x>IO values only in the integration gives results larger by one order of 
magnitude. Because the Coulomb collision cross-section increases as the electron energy 
decreases, using the peak value of electron energy in (73) for a global estimate yields a value that 
underestimates the actual rates of energy losses by collisions with ions. It is therefore clear that 
the value claimed in [2] is unphysical in a dramatic way if the ions and electrons are co-located, 
as indicated by the theoretical model of the EST of f I]. 

It was suggested by EPS that the EST surface should be considered as a strongly coupled system 
with collective behavior, and that particle-particle collisions do not apply. However, this 
suggestion is at odds with the fundamental physics involved. It is true that the electrons in the 
collective state do not “collide” with each other, and if the ions also form a strongly coupled 
plasma, ion-ion collisions are also removed. However, the ions and electrons together do not 
form a strongly coupled plasma! The condition (70) for strong coupling is not satisfied for the 
electron-ion system, because the relative difference in energy (10 keV) between ions and 
electrons is far too large. Therefore, high-energy electron-ion collisions definitely occur. The 
situation can be better visualized with the following analogy. Consider a metallic solid with a 
conduction band or an ionic crystal. The system can be thought of as an ion component in a 
lattice (strongly coupled plasma) with a neutralizing background of electrons. Collisions between 
ions do not occur in this system; only collective excitations (in this case, phonons) are naturally 
present. Yet, if a 10 keV atomic beam is directed at the solid, one is guaranteed to have a lot of 
individual particle collisions with a significant amount of ions ejected from the solid. The case of 
the EST is similar. Ions and electrons would respectively form lattices at near solid density, but 
impacting each other at high velocities. Because the relative kinetic energy between the two 
components is large, particle collisions frequently occur and are correctly described by Equations 
(73) through (81). 

The radiative losses should also be considered. The rate of cyclotron emission can be obtained 7 
from ref. [II] again (p. 68): 

P t =6.21xl(T 17 n„T eaeV) B 2 W/nv (82) 


’ Mlhough solid copper is not a plasma per se, the Rutherford cross-section would still apply at these high energies. 

Strictly speaking, since the electrons are mono-energetic, one could use the radiated power emitted by a single 
electron, using classical formulas such as ( 14.3 1 ) of f 1 3 1. However, this does not affect the discussion. 
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Here, 7 = E e , and this gives a rate of total energy loss by cyclotron emission for the EST: 

Eg* » 6.95 x10 s Br J/s (83) 

If the field is 16,000 T as claimed, the rate of energy loss by cyclotron emission is: 

Eg*- 1.8xl0 14 J/s (84) 

again a difference of 14 orders of magnitude with the rate claimed in [2] (see Table 3). However, 
this would be a gross overestimate of the emitted power, since the highest electron velocities are 
in a region near the outer edge of the plasma shell where the magnetic field is small. To correctly 
evaluate the cyclotron radiative power, the numerical integration scheme was modified to include 
an integration of the cyclotron power. A solution of the dynamical equations was chosen such that 
the maximum field was 16,000 Tesla and the maximum electron energy 10 keV, for the same size 
as the EST configuration given in Table 3. Of course, the corresponding plasma densities and 
total number of electrons were much higher than those derived from the EST mass, following the 
approximate scaling relationships discussed above. The total emitted power for this solution was. 

E'""' = 4xl0 7 J/s (83) 

/;.S7 

while Equation (82), using the maximum values globally, would give 10 22 J/s. Therefore, there is 
a large reduction in emitted power when the correlation between electron energy and magnetic 
field is taken into account. Nevertheless, the solution gives a cyclotron power much larger than 
the one quoted in [2], However, it would become of the same order if the results were scaled to 
match the density obtained from the EST mass. Therefore, the cyclotron power emitted may be 
correct, but the quoted magnetic field value is in serious disagreement with the numerical 
solutions and the scaling laws. 

The total energy can be evaluated as the magnetic field times the total volume of the torus, i.e.. 



Assuming B * 16,000 T, the total energy is 6xlO l? J. Even if the sizes given in Table 3 were the 

actual diameters instead of radii, the energy would be 7xl0 12 J. This is still larger than the value 
mentioned in [2] and listed in Table 3, yet another inconsistency in the claimed results. It can be 
easily verified that the kinetic energy of the electrons is negligible compared to the magnetic 
energy. If the rates of energy losses given in [2] were used, the EST would have a lifetime of 
50,000 years! Using the loss rate from Equation (78), the EST would last about 7000 sec, still an 
extraordinary achievement in terms of plasma stability. However, the true collisional lifetime is 
given by the electron energy lifetime itself, which according to Equation (77), is only on the order 
of a millisecond. 

Therefore, it is clear that collisions between electrons and ions would prevent the EST to achieve 
the very long claimed lifetimes. Of course, this is a problem if the electrons and ions are both 
located in the same regions, which is a key assumption in the model described in [1]. Since this 
model was repeatedly mentioned as an independent description of the physics of the EST, it is 
then reasonable to assume that it is the current state-of-the art in the theoretical description of the 
EST. In fact, as mentioned in Section 4 regarding the ion dynamics, there is no clear description 
of the ion properties anywhere else in the literature [2-3]. The new version of a document [5] 
explicitly mentions the separation of ions and electrons. This is a new addition, made after a 
preliminary version of the present report was made available to EPS. 

Separation of ions and electrons is possible in some cases. As pointed out by Seward [2,3], 
strongly coupled one-component plasmas can organize themselves in shell-like structures, as 
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demonstrated by Gilbert et al [14] for an ion plasma in a Penning trap (i.e., confined by an strong 
external magnetic field). Strongly coupled plasmas are now commonly created in the laboratory, 
although always in the presence of confining field configurations such as Penning and Malmberg 
traps. Strong coupling is indeed a key assumption f of the EST model, thereby implying that the 
ions and electrons form two distinct components, each one at low temperature. In fact, the 
coupling condition T >170 for EST is typical of a solid-like plasma phase, corresponding to 
Coulomb crystal structures [15] rather than fluid-like shells, and it is probably more stringent than 
necessary. Due to the large mass difference between ions and electrons, one could then expect ion 
and electron shells to be separated. However, no one has yet realized such a two-component 
system. Furthermore, the realization of strongly coupled plasma configurations is rather delicate 
and currently requires cooling to very low temperatures in the presence of external confining 
fields. To expect that such structures can be formed from high-energy beams or electric arcs is yet 
another formidable claim that demands extensive proof. Finally, the ion shell would also be 
fundamentally unstable, since the system is globally charged with the relative excess of ions. 
Therefore, the ion shell would not remain separate from the electrons, and again, the collisions 
between electrons and ions would be unavoidable. 

The ion confinement problem remains a key issue with important consequences. For example, 
suppose that the ions reside in a shell interior to the electron shell. If mechanical pressure is 
required to confine the ions, as now suggested by EPS [4], there must be a significant number 
density of neutrals at the edge of the ion shell and therefore within the electron shell. Then, the 
collisions between electrons and neutrals become inevitable, and at this range of electron energies 
(10 keV), the previous classical analysis remains relevant. Therefore, excessive collisional energy 
dissipation (with both ions and neutrals), remains an important problem that prevents any 
possibility of long lifetimes for any EST configurations. 

After mentioning the need for a confining force provided by an external pressure, EPS again 
modified the theoretical model and stated [4] that there is a fixed ion background. The latter 
would be provided by an ionization process of the neutral particles diffusing inside the EST. In 
that case, electron-neutral collisions must be considered. The range data [16] of 10 keV electrons 
in air specifies a figure of 19.7 MeV.cm'/g from which one can derive the electron-neutral elastic 
collision cross-section to be of the order of 10' 20 m 2 . If the neutral density were of the same order 
as the ion density in the EST, the electron lifetime would be 10' 8 sec. If instead, only a partial 
pressure of 1 pTorr is present, the lifetime would be approximately 50 msec. This is still 
significantly less than the 10 6 hours claimed in [2], At this point, it becomes important to evaluate 
the actual neutral partial pressure that would be required to maintain the EST in this scenario. 
This was already done in Section 5, thereby leading to a relationship between the external 
pressure and the magnetic field pressure. Equation (65). Another approach is now used, and 
again, the case of high-energy EST (Table 3) is examined. 

The electrostatic confinement of the electrons (from the excess ionic charge) is responsible for the 
presumed generation of the large magnetic fields inside the EST (16,000 T in this case). From 
energetic considerations, if the field energy is much larger than the kinetic energy of the plasma, 

one should expect a balance between the magnetic field energy density ( B 2 /2(u 0 ) and the energy 

density (or pressure) of the electrostatic field ( e 0 E 2 12). As the former decreases near the edge, 
the latter increases. This indicates that the electric field at the edge of the shell is of the order of 
5xl0 12 V/m. Another evaluation comes from the total charge of the EST, quoted in [2] to be 186 
Coulombs. Assuming a cylindrical shape, the field would be: 


It is not a derived condition, as implied in some of the EST literature. 
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(87) 


(g/2w/? L ) gl . 6x|0 i3 v/m 

2nr b 

These two approximations are in fairly close agreement. One can then consider the steady-state 
dynamics of the ions near the edge of the shell. The magnetic field vanishes in that region, and 
the electric field provides a large acceleration, while the collisions with the neutrals provide a 
frictional force, i.e.: 

dv eE 
dx mj 

where * is the distance along the ion path. Since the collision frequency between ions and 
neutrals also depends on the ion velocity (v,„ = N n a, n v), an approximate solution can be found 


of the form: 


V 

2 


— a(|-c-* /A ) 

ntj 


(89) 


where A = (IN n c ,„) 1 is the mean free-path. Note that as the magnetic field vanishes towards the 
edge of the shell, the Larmor radius becomes large, and the collision mean free-path becomes the 
relevant length scale'. The cross-section is typically of the order of 10“ nr, and if the gas 
pressure is lpTorr at 300 K, the mean free path is 1,500 meters. The ions would therefore stream 
away from the EST at extremely high velocities without stopping, thereby invalidating the 
concept of partial vacuum confinement. 

If the gas density is that of the presumed ion background inside the EST, the mean free path 
becomes 300 pm, which is a small fraction of the shell thickness. It would appear then that this is 
the range of neutral density required. A more accurate determination is shown below. It is 
interesting to also point out that the velocity gain during this distance is (according to (89)), about 
7x10* m/s! This value is approximately the same as the one obtained from the ExB drift 
velocity, using peak values of both fields. 


It was previously mentioned how energetic considerations dictate that the neutral gas pressure, in 
order to confine the ions, should be of the same order as the electrostatic energy density in order 
to confine the ions; i.e., of the order of 2.5xl0 9 atmospheres for the EST described in [2] and 
Table 3. A similar result can now be obtained from kinetic considerations. Equilibrium is 
achieved when the momentum fluxes at the EST interface can be achieved. The ions inside the 
EST are accelerated to high velocity from the electric field, and the velocity at the interface can 
be obtained from (89) using an average over a mean free path. The representative distance is 
therefore A / e = A / 3 , and the average ion velocity at the EST boundary is: 


The ion momentum flux across the interface is Nj(mjV)v and must be balanced 
corresponding flux from the neutral background (see Figure 15), i.e.. 

N, (m i v)v = N n (m n v 2 x ) = N„kT„ 


(90) 
by the 


(91) 


where the statistical average ( v 2 ) over the neutral particle distribution has been used. Using (87), 
the required neutral pressure then becomes 9xlO' 3 Pa, or 900 million atmospheres; again of the 
same order as the magnetic field pressure. 


Even using the Larmor radius for peak magnetic field does not actually alter the conclusions. 
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Figure 15: Schematic of kinetics at the EST interface. 


More information can be gained from this elementary analysis. If the scenario of confinement by 
external pressure is abandoned and one relies instead on the ionization approach to maintain a 
fixed ion background, the actual rate of ionization can now be estimated. Again, using (90), the 

average ion velocity at the interface can be evaluated at = 5xl0 8 m/s. Using the actual 
dimensions of the EST (Table 3) and the inferred density, one can then evaluate the rate of 
particle loss. This is found to take the staggering value of: I0 32 particles/sec. Since it takes at least 
13.5 eV to ionize the gas, this corresponds to an energy loss rate of approximately: 

E'est ~ 2.3 xlO 14 J/s (92) 

again a difference ol 14 orders of magnitude from the claim made in [2], Since this energy loss 
directly affects the kinetic energy of the electron component (=10 21 electrons x 10 keV), it is 
found to correspond to an EST lifetime of less than I0 8 sec, compared to the 10 6 hours lifetime 
claimed in [2], 


In a more recent document [8], Chen claims that the EST confinement is not incompatible with 
the Virial theorem. However, this statement deserves to be investigated in more detail, since the 
conditions necessary for this to occur still make the concept of long-term stability of a free EST 
nonfeasible. One must first start with the complete description of the Virial theorem: 

^~r = 2 ( E kiii +E thr) + Eflec +E maii ~ jd* X X a ^— - j dS p X (f ( + 7 '°^ ) ( 93 ) 

where / is the moment of inertia of the plasma, E kin is its kinetic energy, E, hr is the thermal (or 

random kinetic) energy, and E elec ,E inaK are the energies of the electrostatic and magnetic fields 

respectively. If the plasma is a collection of particles of various types a , the following 
definitions apply: 

I = ^jm a n a r 2 d*x (94) 

o y 


E k „ , = Y,jm a n a 
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d, I 


E thr = X J m <? n O ((V - ft ) 2 y*x 
a V 

where u = dx / dt is the mean plasma velocity. The field energies are of course given by: 

E e ,er=\ 'd*X$£- 

v z 


(95) 
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The radiation momentum density vector G and Maxwell stress tensor are given by: 


G = £ () ExB 
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For a sufficiently collisional plasma, the mechanical pressure tensor p a> 3 only has a diagonal 


component p = nkT . For a plasma in steady state, the moment of inertia is zero, i.e., the plasma 

neither expands (d 2 l /dr >0) nor contracts ( cl 2 / / dt 2 < 0). Assuming also that the plasma is at 
rest ( E kin = 0) and using the relationship between pressure and thermal energy, the condition for 


stability becomes: 

3 pV + E e , ei . + E max =jd'xx (X ^ + \dS p x a T al} +4 nR 2 p exl (101) 

where the boundary term for the mechanical pressure was computed assuming the plasma volume 
as a sphere of radius R . The first term on the right-hand-side of (101) describes the effect of the 
radiation of electromagnetic energy: it physically corresponds to the “recoil” from the emission of 
photons. The second term describes the electromagnetic stress provided at the boundary by an 
external (confining) electromagnetic field. The third term is the effect of a mechanical pressure. If 
the plasma does not radiate and if there are no external fields or pressure, the RHS of (101) is 
identically zero. Since all terms on the LHS are positive definite, the stability condition cannot be 
satisfied, and the plasma is unstable (it expands). This is the origin of the impossibility of self- 
confinement of a plasma. Chen claims that these RHS terms are not necessarily vanishing, and 
can therefore provide a method of confinement. However, as we shall see, these terms would 
need to be very large. 


Let us first consider confinement by radiation. According to [1], the Poynting vector due to 
coherent perturbations of the field at the surface of the EST is: 

S = UE 0 8B e e 2mw)r (102) 

It is of second-order with respect to the perturbations and therefore very small. Furthermore, for 
stable configurations, Im(co)<0; and the radiation decays exponentially in time. The same can be 
said for the momentum density vector G . Chen uses the same coherent perturbations for the 
contribution of the Maxwell stress tensor at the boundary, and these are also of the same type. 
Therefore it is difficult to see how second-order perturbations, decaying exponentially in time, 
would contribute to long-term stability of the EST! Generally speaking, there is a fundamental 
problem with invoking radiative processes to provide a confinement mechanism. Physically 
speaking, conservation of momentum during the emission of radiation provides the mechanism 
for preventing the plasma particles to escape the system. Thus, if there is a net radiative 
momentum outward flux, the recoil of the plasma keeps it momentarily stable. However, the 
momentum of a photon is a rather small quantity ( hv / c ). In order to keep the EST confined, a 
tremendous amount of power would need to be radiated (c times the momentum flux!). The 
radiative losses would therefore lead to an extremely rapid collapse of the plasma. 
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The other term mentioned by Chen is the contribution of the external mechanical pressure. 
However, it is easy to see that this pressure must be large. Neglecting the radiative terms, one can 
see that in order to have confinement: 


Pex, = P + ~ 


,E 2 


B 2 


(103) 
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Therefore, the external pressure must be of the same order as the energy density stored inside the 
EST, as discussed previously. It is difficult to imagine how Chen proposes to confine a high- 
energy EST with an internal magnetic pressure of 1 billion atmospheres in partial vacuum. 


It is clear that none of the documented or potential explanations provided by EPS for the claims 
of stable energy storage is physically possible. So far, a valid theoretical model that can account 
for the experimental observations remains elusive. The extraordinary claims regarding the 
properties of the EST as well as the theoretical procedures utilized by EPAS are at odds with the 
known relevant physics. One could be tempted to assume that, although the theoretical models 
provided by EPS are in error, there could be experimental evidence of some structure. However, 
based upon the additional documentation provided by EPS, there are also some problems with 
this type of evidence, as described in the next section. 


8. Experimental Observations 

Among the additional documentation furnished by EPS was a report to the defense Special 
Weapons Agency (DSWA) [6], which included a slightly more elaborate description of the 
experimental procedure used to create ESTs. In the 1999 version of an EPS document [3], EST 
initiation is described in terms of electron beam injection, while in the 2000 version [5] an 
electric arc is used. There is also a mention in [6] of generating the ESTs by using an exploding 
wire. The electric arc approach is described in [6]. In a tabletop experiment, a high-current arc 
discharge (400 amps) between two carbon electrodes is initiated in a low-pressure (0.001 1 atm) 
Nitrogen atmosphere. The experiment lasts 10 seconds and is characterized by electrode material 
vaporization and the melting of parts of the structure. Pictures of the experimental set-up show 
significant thermal damage. The “proof’ of the presence of ESTs consists of a picture of a small, 
luminous ring-like object seen at some distance from the arc. It is not clear whether this is an EST 
or another more benign structure, since there is no measurement of some of the key 
characteristics of the EST.. For example, plasmoids (neutral plasma structures with toroidal 
configurations), have often been generated [17] and have demonstrated (relatively) good stability 
properties, although certainly not on the level claimed in [2-6]. There is also repeated mention of 
the charge neutrality of the EST, since it is not affected by the electric field. This would then also 
imply that there is no significant energy storage or magnetic field trapped within the torus, since 
the two are closely related. 

An energy loss rate of 1.89xl(T 7 joules/m 2 of surface area (assuming that the authors meant 
watts/m 2 ) is quoted [6] and claimed to be an important result that demonstrates the stability of the 
EST. It is worth looking at this in more detail. The luminous power of the EST is obtained in [6] 
by visually comparing it with that of a Light Emitting Diode (LED) and is (taking into account a 
geometric factor) of the order of 0.0026 watts 7 . The EST has a major radius of 3.5 mm and a 

minor radius of 1 mm, thereby giving it a shell surface area 1.38X10 -4 m 2 . Dividing the total 
emitted power by the surface area, one should then obtain an intensity of 19 watts/m 2 . The result 

of 1.89x10 7 joules/m' claimed in [6] is obtained by multiplying (instead of dividing) by the 


Again, the units used in [6J were Joules for a rate, instead of watts, or J/s. 
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surface area! This “important result” of low energy loss rate is therefore incorrect by 8 orders of 
magnitude. This error has since been corrected by EPS. 

In the same reference ([6], p. 1 1), a “demonstration” EST is described with characteristics similar 
to the one “observed” in the arc experiment, at least as far as the electron energy (1 10 eV, of the 
order of the energy acquired in the cathode layer) and the ambient atmosphere are concerned. 
However, the major radius is now 10 cm, while the radiative power is now 0.00046 J/s. It is not 
clear how a total radiative power 10 times smaller than the experimental one could be obtained 
when the surface area is 3000 times larger than the (assumed) experimental EST. 

It is then further stated that the number of ions is approximately 1% greater than the number of 
electrons. This is in disagreement with the values quoted in the same document for the analysis of 
the experiment. In this analysis, Chen quotes a value for the relative charge excess: 

/-l = 7xl(T 5 , which is significantly different from the 1% mentioned elsewhere. It is also 

worth looking at this experimental analysis in more detail. A key assumption in Chen s analysis is 
that the background ion density is equal to the ambient density, which is approximately 

n * 2.5 x 10 22 m \ Using this value for the electron density and the given EST dimensions, this 
corresponds to a = 890, which is close to the value (a = 1100) quoted by Chen. Another 
assumption is that the maximum electron energy is 200 eV. Using this value and the computed 

value of a with the scaling relation (37), one finds that /-l = 8.7xKT 5 , again close to the 
value quoted by Chen. However, there are two other important pieces of information concerning 
the experiment that are important to mention: 

. The observed ESTs appear to be charge-neutral, since they do not appear to be affected by the 
electric field near the electrodes. 

• The EST’s are affected by gravity, and are seen to follow parabolic trajectories. 

The latter can provide some information about the typical acceleration to which the EST is 
subjected. Using a characteristic dimension of L = 10 cm (the average interelectrode separation) 
and an EST lifetime of 0.6 sec, the acceleration must be at most of the order of: 

%,v4 x0 ' 5m/s2 (104) 

t~ 

Note that this corresponds to an average mass density in the EST only 5% larger than the ambient 
density, which is consistent with Chen’s assumption. However, the acceleration in an electric 
field is simply given by: 

a e , «.-(/-!)—£ ( 105 ) 

m N 2 

where m N ^ is the molecular mass of nitrogen. With the applied voltage of 200 Volts used by 
Chen and an interelectrode spacing of 10 cm, this leads to: 

a elec ~ 5 x 1 O' 5 ’ Hl/s 2 ( 1 06) 

In order to match the observation that the EST is “charge-neutral” (i.e., is not affected by the 
electric field), the acceleration from the field should be much less than the one due to gravity. 
However, according to the parameters deduced by Chen, it is 6 orders of magnitude larger. 

Therefore, it is clear that this document provides no reasonable evidence supporting the existence 
of the high-energy EST (i.e., one with a significant electric charge). However, a completely 
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charge neutral structure (and therefore probably not an EST) is not ruled out. It is also interesting 
to point out that a rate of energy losses by collisions was given for the very high energy EST 
described in [2], which is extremely low to support lifetimes of the order of several years. It is not 
clear how this value was obtained. As demonstrated earlier, if the same theoretical model [I] used 
to describe the properties of this EST is utilized, the collisional losses mentioned in [2] are 18 
orders of magnitude below the expected levels. If the ions and electrons are separated as claimed 
in the more recent documentation [5], there are no collisional losses; however, the ion 
confinement remains a daunting problem. If neutrals must exert a confining pressure on the ion 
shell inside the torus, then electron-neutral collisions are unavoidable and the energy losses are 
again large. If the losses are estimated from the experimental data (assuming that the observed 
structures are ESTs) and using the rate of energy loss per unit of EST surface area the results are 
underestimated by at least 8 orders of magnitude as shown above. Despite claims to the contrary, 
there is no justification for assuming a low rate of energy loss of the EST. 

In a recent meeting between MSE, EPS and MIT [7], additional experimental evidence was 
presented. Primarily, this consisted of photographs of toroids produced in arc experiments and in 
exploding wire experiments. The latter method appeared to produce a much larger number of 
toroids (the assumed EST) than the arc method. An interesting parametric study was also 
performed by EPS with respect to the appearance of the toroids. In that study, a photographic 
image was simulated, assuming a toroidal plasma shape and various rates of rotation along major 
and minor axes. This study was then able to provide the range of rotation speeds required for the 
photographic image to always appear as it does in the experiments. This study was of great 
interest, because it provides the first example of valuable data against which various physical 
models of the plasma toroids can be compared. The study is definitely a step in the right 
direction. One key difficulty in the experimental investigations of the EST is the lack of 
diagnostics. This is of course related to the lack of reproducibility. Since the production of ESTs 
is an inherent random phenomenon, it is difficult to predict the time and location at which the 
EST is formed, and where advanced diagnostics could be applied. Therefore, it appears critical to 
proceed in two directions: 1) systematically investigate the conditions for maximal rate of EST 
production; and 2) perform systematic statistical analyses of the produced EST. The difficulties 
associated with EST characterization are duly appreciated. Nevertheless, it is strongly urged that 
further progress in this direction be undertaken, in order to provide useful experimental data to 
test all possible models of the observed phenomena. 


9. Conclusions 

Based upon the analysis of all the provided information, a number of conclusions can be made 
regarding the EST concept as proposed and observed by EPS. It is important to distinguish 
between the theoretical models and the experimental observations. It is undeniable that EPS does 
observe luminous toroidal structures in its experiments. However, the nature of these structures is 
far from clear. For example, one could simply assume that they are luminous fluid vortices 
(“smoke” rings) containing small particles at high temperature (the result of electrode or wire 
vaporization). Photographic snapshots could not provide enough information to eliminate this 
tentative explanation, but the study on rotation rates indicates that this is extremely unlikely (thus 
the usefulness of such studies, as mentioned earlier). However, a rigid structure cannot be ruled 
out. Other various complex, exotic and likely models can be proposed, as well. The cold-fluid, 
multishell model with a high magnetic field stored within the toroidal shell as proposed by EPS is 
rather extreme and (as described in this report) highly unlikely. However, there may be some 
hope of finding the correct physical model. Once this is done, potential applications can be 
accurately assessed. 
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MSE’s conclusions regarding the analytical and experimental data published by EPS and MIT 

regarding the EST concept are as follows: 

( 1 ) The cold-fluid equations for the electron gas and their numerical solutions (described in [ I ]) 
have been verified. These equations are applied to a toroidal plasma shell in the 
approximation of a very large torus diameter ( R T » r h ). The electron fluid rotates to create 
a magnetic field inside the torus, while being confined by the electrostatic field created by a 
excess of positive charges in the plasma. MSE’s numerical solutions agree with Chen’s on 
this point. 

(2) However, the model in [1] does not describe the complete dynamics of the system. Within 
the context of that model, it is found that the ionic component is unstable in the presence of 
the required electrostatic field. Therefore, one cannot at the same time realize a 
configuration that confines both electrons and ions without invoking an external 
confinement mechanism. An external confining field would need to be of the same order as 
the field being stored inside the torus, thereby limiting practical applications of the concept. 

(3) The discrete model proposed by EPS contains various fundamental errors. In the limit of a 
large number of electrons, the correct equations were provided by MSE and are consistent 
with the cold fluid model. Again, the ions cannot be naturally confined. 

(4) An external mechanical pressure is required to confine the ions. This was mentioned by 
Chen after MSE provided a preliminary copy of the present analysis to EPS, but a small 
(“partial” vacuum) pressure was deemed necessary [4]. However, it can be seen from simple 
arguments and from the classical Virial theorem, that the confining pressure would need to 
be at least of the same order of magnitude as the energy density contained in the EST. 
Therefore, there is no possibility of confinement of high-energy ESTs in a partial vacuum. 

(5) The other contributions to the Virial suggested by Chen [8], namely radiative emission and 
perturbations of the electromagnetic stress tensor at the boundary due to coherent 
fluctuations are not relevant. Physically speaking, it is pointless to use the radiative emission 
to prevent plasma expansion, since this would imply a very high rate of energy loss, thereby 
leading to a rapid collapse of the structure. 

(6) The need for an external pressure (mechanical or electromagnetic) of the same order as the 
contained energy density is critical to the assessment of feasibility and usefulness of the EST 
concept. Any potential for large energy storage capability with no external field is negated, 
since the requirements for an external magnetic field are replaced by an equally impractical 
requirement for large neutral pressure. 

(7) Using the current theoretical model of the EST, it has been shown that the mixture of ions 
and electrons would be the subject of excessive losses through Coulomb collisions and 
possibly cyclotron radiation. This energy loss rate is sufficient to prevent “long-term” 
stability of the plasma, even in the presence of an externally applied confining field. A 
modification of this model using separation of ions and electrons into separate “shells” 
would prevent these collisional losses but is not possible without a plausible mechanism for 
ion confinement. Since an external pressure would be required to confine the ions, electron- 
neutral collisions then become unavoidable and lead to similarly large and debilitating rates 
of energy losses. 

(8) The mention by EPS [4] of suppression of electron-ion collisions due to the strong coupling 
is not valid. Since the relative velocity between electrons and ions is large, the electrons and 
ions together do not form a strongly coupled plasma. 


(9) After further discussion [4], EPS stated that the fixed ion background was real and was the 
result of a continuous ionization process of the background neutrals. This process was also 
examined and found to be fundamentally Hawed, leading to even greater rates of energy 
losses that are incompatible with any claim of long term stability. 

( 10) The numbers presented in the literature [2-6] are often inconsistent and at times unphysical. 
Furthermore, the analysis of the experiment presented in [6] is inconsistent with the 
observations. There is also no measurement of the essential physical properties of the EST, 
so other interpretations of the experiment cannot be eliminated. 

Based upon these results, one must conclude that the EST concept (as currently proposed by EPS 
and MIT) is fundamentally flawed. None of the claims of extraordinary energy storage capability 
and absolute stability can be scientifically substantiated. Although the present analysis used the 
example of a high-energy EST for demonstration purposes, the conclusions affect EST 
configurations on all scales. Various tentative models for ion confinement and EST stability have 
been suggested by EPS, but none of these satisfies basic physical constraints. These claims were 
made before the problem of ion confinement had been seriously addressed. It is absolutely 
essential that this problem be resolved before long-term stability and energy storage properties 
can be accurately evaluated. 

Although the documentation concerning the experimental evidence is limited, it also appears that 
a thorough analysis of the experimental data has yet to be performed. The experiment described 
in the provided reference [6] utilized inadequate diagnostic procedures. However, it should be 
pointed out that this is not necessarily the responsibility of the investigators, but rather the result 
of the chaotic nature of the phenomena and possibly the lack of adequate funding. This analysis is 
the most complete to date and used the very latest data and models suggested by EPS. It also 
exposed a number of fundamental problems, which had not been previously identified. 

The principal objections to the EST concept as currently proposed by EPS concern the theoretical 
model, as well as claims of very large energy storage and stability properties as made by its 
investigators. Unfortunately, these claims divert attention from other possible physical 
explanations tor the experimental observations. As emphasized earlier, the observed production 
by EPS of plasma structures in exploding wire and arc experiments is undeniable. Although this 
is not the first time that such plasma structures have been produced [18], the EPS observations are 
relatively new and can potentially benefit from better diagnostic procedures. The study of “Ball 
Lightning” [19] phenomena contains many similar observations, along with a large variety of 
theoretical models. Some of these models could also be invoked as an explanation for the 
experimental observations made by EPS. At any rate, the theoretical work done to date by the 
EPS/MIT team lacks completeness and accuracy. For example, the ion confinement problem was 
never mentioned before discussion with MSE, yet it is critical to long-term stability. Various 
mechanisms were proposed to address the stability problem, but none was found to be physically 
feasible. The current direction in theoretical development by MIT now seems to include a two- 
fluid model, yet it is not clear what could be expected from this new development. Instead, it 
seems clear from a look at the Virial Equation [e.g., Equation (101)] that the real problem 
concerns the lack of a negative energy term. For example, this is supplied in astrophysics by 
gravity such as in the cases of self-gravitating plasma rings [20]. In a degenerate plasma, the 
exchange energy plays a similar role, and this allows an ionic crystal to be stable. This would be 
an alternative possible approach to the classical fluid model. Of course, such a system could not 
allow the storage ot energy in a magnetic field, since both the electrons and ions together would 
need to be strongly coupled (i.e.. no relative velocity). This model could have some applications, 
although it would not be an EST (i.e., no magnetic field). 
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It should now be clear that there is no physically feasible solution for the EST concept which by 
definition includes high loop current, non-neutral plasma, and trapped magnetic field. As 
proposed, the energy storage applications are not possible. However, the possibility of charge 
neutral, metastable structures (lifetimes of the order of several milliseconds) is intriguing and can 
have a number of applications. Chen proposed to use such structures for propulsion, since their 
long lifetime ( I ms) might allow them to be accelerated to high speed. With the structure being 
charge-neutral, the only way to apply an accelerating force (as proposed by Chen) would require 
a gradient of magnetic pressure. However, to generate a strong gradient across the EST dimension 
(< 1cm) with an applied magnetic field would be a serious challenge. 

MSE believes that the pursuit of stable (or rather, “metastable”) and energetic plasma structures is 
a worthy and very important endeavor. Research in this direction should be pursued, but with 
much caution and under strict peer-review. A more determined and focused research program on 
theoretical model development should be a high-priority item. Despite the many problems with 
this EST concept, it is recommended that independent theoretical (using different models) and 
experimental investigations of the feasibility of using various plasma structures for propulsion or 
energy production should be initiated by NASA. 
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